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Abstract 

The goals of this paper are to show the following. First, Grover's 
algorithm can be viewed as a digital approximation to the analog quan- 
tum algorithm proposed in "An Analog Analogue of a Digital Quantum 
Computation", by E. Farhi and S. Gutmann, Phys. Rev. A 57, 2403-2406 
(1998), quant-ph/9612026. We will call the above analog algorithm the 
Grover-Farhi-Gutmann or GFG algorithm. Second, the propagator of the 
GFG algorithm can be written as a sum-over-paths formula and given 
a sum-over-path interpretation, i.e., a Feynman path sum/integral. We 
will use nonstandard analysis to do this. Third, in the semi-classical limit 
ft — > 0, both the Grover and the GFG algorithms (viewed in the setting 
of the approximation in this paper) must run instantaneously. Finally, 
we will end the paper with an open question. In "Semiclassical Shor's 
Algorithm", by P. Giorda, et al, Phys. Rev. A 70, 032303 (2004), quant- 
ph/0303037, the authors proposed building semi-classical quantum com- 
puters to run Shor's algorithm because the success probability of Shor's 
algorithm does not change much in the semi-classical limit. We ask the 
open questions: In the semi-classical limit, does Shor's algorithm have to 
run instantaneously? 



1 Introduction 

This paper attempts to answer the question: what is a quantum algorithm? In 
particular, what is Grover's algorithm? The views which we put forth in this 
paper are that quantum mechanics is an analog, continuous-time theory, and 
Grover's algorithm is a digitization of analog quantum mechanics. In particular, 
it is a digitization of the analog quantum search algorithm proposed by [7]. 
The digitization is performed in such a way so that the success probability of 
the analog search algorithm is well preserved. We are not certain that this 
perspective of Grover's algorithm is the correct one, but it is certainly true that 
the theory of quantum mechanics describes quantum evolution in continuous 
time. 

The main goals of this paper is to show some relationships between the 
Grover algorithm and the analog quantum algorithm and elaborate more on the 
analog algorithm. We will call the analog algorithm proposed by [7] the GFG 
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(Grovcr-Farhi-Gutmann) algorithm. The three things which we will show in this 
paper are the following. We will show that Grover's algorithm can be viewed 
as a digital approximation to the GFG algorithm. This interpretation comes 
about via the Lie- Trotter product formula. Second, We will use the Lie- Trotter 
product formula and derive a sum-over-paths formula for the propagator of the 
GFG algorithm. We will then give the sum-over-paths formula a Feynman path 
summation (integral) interpretation. This requires some nonstandard analysis. 
Finally We will investigate the semi-classical limit behaviors of the GFG and 
the Grover algorithms (in the digital approximation setting of this paper). This 
will lead us to ask an open question for Shor's algorithm (for Shor's algorithm, 
sex: [19] and references within). We will assume that the reader is familiar with 
both the Grover and the GFG algorithms, Feynman path integrals, and a small 
amount of nonstandard analysis. For the readers who do not fall into the above 
category, we will lightly sketch the details when we go into those subjects. 

In a few sentences, the GFG algorithm runs as follows. Given a Hamiltonian 
H w = E\w){w\, drive the Hamiltonian H w into H = H w + H s where 11 H s = 
-^(l s )( s l — -0- Evolve the initial state \s) via the evolution e ~ ltH / h for some 
time t = t , then \{w\e~ lt H / h \s)\ 2 1, i.e., the algorithm finds \w). 

It turns out that it is possible to derive an interesting relationship between 
the GFG algorithm just described and Grover's algorithm. This relationship 
comes from the Lie- Trotter product formula and it naturally leads to the inter- 
pretation that Grover's algorithm is a digital approximation of the GFG algo- 
rithm. Using the Lie- Trotter product formula, we can break the above evolution 
into 



By setting t to the appropriate value, the second line of expression 1.1 becomes 
a fc-product of the "Grover search engine" 1,2 . In general, setting t to the appro- 
priate value destroys the Lie- Trotter approximation in 1.1 in the sense that the 
error is no longer bounded by inverse powers of k. The situation is remedied by 
setting k to the proper value. Doing so, the success probability of finding \w) 
in the GFG algorithm is well preserved (with error bounded by inverse powers 
of k) in the approximation. The analysis of the former error was previously 
done in [21]. We can attempt to derive a satisfying (1) bound on the former 
error, see [21] (we will briefly outline this in section 4). Doing this leads to the 
results in [21], but this is not our goal. Our goal in this paper is to concentrate 
on the analysis of the latter error. We wish to show that as the size of the 
database to be search approaches infinity, the success probability of Grover's 
algorithm approaches the success probability of the GFG algorithm, i.e., the 
error between the two probabilities is bounded by inverse powers of the size of 

11 In [7], H s = E\s)(s\. Hence, our evolution differs from the one in [7] by a phase. 
1,2 This term was coined by [4], see equation 2.4. 
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the database. Because of this, we take on the view that Grover's algorithm is a 
digital approximation of the GFG algorithm. In this sense, we can view Grover's 
algorithm as a special case of the Lie- Trotter product cut-off approximation for 
the GFG algorithm. By special case, we mean that the approximation is good 
for the success probabilities with input vector \s) and output vector \w), i.e., 
the approximation is not global and it does not preserve the wave functions. We 
will show this section 4. Given the current state of the art, whether the Grover 
algorithm is an approximation of the GFG algorithm or the GFG algorithm is 
an approximation of the Grover algorithm is a matter of taste. We take on the 
former interpretation because quantum mechanics is an analog theory and the 
second line of 1.1 becomes Grover's algorithm even though it only approximates 
the success probability of the GFG evolution. 

Those who are familiar with Feynman path integral techniques should rec- 
ognize that this is an opportunity for a sum-over-paths interpretation. For 
more on Feynman path integrals, see [1, 2, 3, 8, 9, 13, 22]. The traditional 
non-relativistic quantum mechanics interpretation of the propagator K (x, y, t) 
is that it is the probability amplitude of a particle that starts at position x at 
time zero and ends up at position y at time t. The Feynman path integral inter- 
pretation of the propagator is that for any classical path that starts at x at time 
zero and ends at y at time t, there is an amplitude associated with that path 
and the total amplitude K (x, y, t) is the sum of the amplitudes of all paths. 
With physicists' rigor, the Feynman path integral can be derived by breaking 
up the evolution operator via the Lie- Trotter product formula, interpreting the 
limit as time slicing the physical evolution and then going to the continuous 
time limit. From a mathematician's point of view, this process must be done in 
imaginary time since in real time, there is no measure for which the Feynman 
path integral converges (see [2, 3]). 

Extending this into quantum computing, the propagator K (j, k, t) would 
be the probability amplitude of starting at state j at time zero and ending up 
at state k at time t. A traditional Feynman path integral interpretation for 
this propagator is somewhat tricky because the set of "positions" for j and k 
is discrete and finite. A more serious problem for the propagator of the GFG 
algorithm is that when time sliced, the summand contains the term 5j m ,( m+1 , 
which is 1 when the state of the system is \l m ) at the m th time slice and is \l m +i) 
at the to + 1 st time slice, and zero otherwise. This term is meaningless in the 
continuous time limit. It is for this reason that we use nonstandard analysis to 
give the GFG propagator a sum-over-paths interpretation. Using nonstandard 
analysis on Feynman path integrals is not a new concept, see [1, 14, 15, 16, 17, 
18]. For more on nonstandard analysis, see [1, 5, 12, 13, 23]. The amount of 
nonstandard analysis that we will use is minimal. We will use the nonstandard 
reformulation of limits. Given a sequence r„ such that 

lim r n = r, (1.2) 

the nonstandard equivalent of equation 1.2 is the following. Given any infinite 
integer w€'N, where *N is the nonstandard extension of the natural numbers 
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N, r w is infinitcsimally close to r, and st(r^) = r. The operation st known as 
the standard part. 

Using nonstandard analysis, we will formulate a sum-over-paths interpreta- 
tion for the propagator of the GHG algorithm as follows. Fix an infinite w£*N, 
time slice [0, t] into the set of discrete 13 times 

T = {0,t/u,2t/iv,...,t}. (1.3) 

For each path 14 P : T — > {1,2,..., N} 7 where N is the size of the data base 
to be searched in Grover's algorithm, we associate an amplitude to that path. 
The propagator is then infinitcsimally close to the sum 1,5 of the amplitudes of 
all such possible paths. Further, this has to be true for all w G *N. We will 
show this in section 5. 

Finally, we will investigate the semi-classical behaviors of the Grover algo- 
rithm and the GFG algorithm in the setting of this paper. We will show that 
in the semi-classical limit, both algorithms must run instantaneously. We must 
point out that if one day we figure out how to build a quantum computer to run 
Grover's algorithm, we do not know whether Grover's algorithm will behave this 
way in the semi-classical limit. This is because in the setting of this paper, we 
think of the Grover algorithm as a digital approximation of the GFG algorithm, 
which may or may not be the appropriate point of view. Recently, the authors in 
[10] proposed building semi-classical quantum computers for Shor's algorithm. 
They showed that in the semi-classical limit, the success probability of Shor's 
algorithm does not change much. At this point, we ask an open question. In 
the semi-classical limit, does Shor's algorithm have to run instantaneously? 



2 The Grover Algorithm 

In this section, we will outline Grover's algorithm. For more details on the 
algorithm, see [4, 7, 11, 19]. The material in this section is taken from [4]. We 
are given a function / : {0, . . . , N — 1} — > {0, 1} such that there exist a w with 
f(w) — 1 and f(a) = for a ^ w. The goal is to use a quantum computer to 
find w. 

Grover's algorithm is as follows. Let Uf be the unitary operator defined by 
U f \a) = (-l) /(a) |a> = I - 2P W , P w = \w)(w\, (2.1) 

and let 

Define U s by 2 - 1 . 

U s = 2\s)(s\ -I = 2P s -I, P a = \s)(s\, (2.3) 

1,3 Hyper-discrete to be more precise 

1 - 4 Internal path to be precise 
1 ' 5 Internal sum to be precise 

2 - 1 In [4], U s is denned as U s = I - 2\s)(s\. 



N-l 



\= E i°>- 



a=0 



(2.2) 
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and let 



U G = U s Uf 



be the "Grover search engine" 2 2 . Let 



|r> 



1 



->-^l 

} y/Nj' 



(2.4) 
(2.5) 



and let V be the subspace spanned by {\w), \r)}. 

Theorem 2.1. The subspace V is an invariant two-dimensional subspace of 
Ug- Further, with respect to the orthonormal basis {\w), \r)} in the invariant 
subspace V, Uq admits the unitary matrix representation 



U 



N-2 
N 



2^fN=l 
N 



N 
N-2 
N 



cos 9 — sin b 
sin 9 cos 9 



sin 



Proof. See [4]. 

Theorem 2.2. LetU G ,\s), 

k 

U k \s) 



(2.6) 
□ 



be as previously defined, then 



cos f — sin b 
sin 9 cos 9 



I N-l 



cos (k9 + a) 
sin (k9 + a) 



a = cos 



1 



The probability of reaching the state \w) after k iterations is 
? fe = \(w\U k \s)\ 2 = cos 2 (k9 + a). 



(2.7) 
(2.8) 

Proof. See [4]. □ 

Basically, theorems 2.1 and 2.2 say that the "Grover search engine" rotates 
the vector \s) in the subspace V to the vector \w) after k iterations. 



Further, at k = ^-,7 k w 1. 



3 The GFG Algorithm 

In this section, we will outline the GFG algorithm. For more details on the 
materials in this section, see [4, 7]. The GFG algorithm is similar to the Grover 
algorithm in the sense that they are both search algorithms. As stated in [7], 
the GHG algorithm solves the following problem. Given a Hamiltonian H w = 
E\w)(w\ where \w) is a basis vector (the Hilbert space has dimension N), find 
\w). 

The GFG algorithm solves the above problem as follows. Drive the Hamil- 
tonian H w into the Hamiltonian 31 , 

H = H W + H s , (3.1) 

2 2 The term was coined by [4] 

31 For our purpose, we will want H = H w + H s = H w + H s — EI. Hence the notation H 
and H s . 
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where 



Let 



JV-l 



a=0 



H s = E\s)(s\, |*> = JL £ |a>. 
1 

|r) = 



w) 
N 



(3.2) 



(3.3) 



then, as in the Grover algorithm, the two dimensional subspace V spanned by 
the orthonormal basis {\r), \w)} is an invariant subspace of the evolution under 
H. In other words 3 2 , 



-itH 



(3.4) 



In the subspace V with respect to this basis, the evolution admits the matrix 
representation 



itH/h 



= e - 



-y/l^ian^ cos^+iysin^ 



(3.5) 



For t > and with the initial state \s), the system evolves as 



4>(t) 



itH/h 



Eyt . . Eyt 

y cos — i sin — — 

h h 



Eyt, 



\w) + \J\-y 1 cos -^-|r) 



(3.6) 



The probability of measuring \w) is then 

m = \(«>w)\ 2 



He -^/ ft|s)|2=sin 2^ +y 2 cos 2^* ; y 



• (3-7) 



Hence, running the system for time t = will measure \w) with probability 

1. 

For our purposes, we want to work with the Hamiltonian 

H = H W + H S -EI = H W + H S = E (\w)(w\ + \s)(s\ — I) = E (P w + P s - I) 

(3.8) 

The evolution of the Hamiltonian H differs from the evolution of H by a factor 
of e lEt / h . 

Theorem 3.1. Let 

H = H W + H S -EI = H W + H S = E (\w)(w\ + \s)(s\ - I) = E (P w + P s — I) . 

(3-9) 

The space V is an invariant subspace of the evolution e ltH / h . In the orthonor- 
mal basis {\w), \r)}, e - ltH / h admits the unitary matrix representation 



-itH/h _ 



-y/l^isin^ cos 2fi + iy sin ^ 



,y : 



(3.10) 



3 - 2 In the quantum computing literature, the value of ft is usually taken to be 1. Because we 
will be taking semi-classical limits, we will not take ft = 1. See section 6. 
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For t > and with the initial state \s), the system evolves as 
<j){t) = e- itH ' h \s) 



Eyt . Eyt 
y cos — i sin --— 

h h 



\w) + \J\ - y 2 cos j • 



(3.11) 



The probability of measuring \w) is 



0>(t) = |H<^)| 2 = K^le-^ls)!- - sin- — + y-cos- — , y = 



,2 • 2 ^ , 2 2 ^ 

■ - — + rcos z ^, 



(3.12) 

Hence, running the system for time t = 7Th ^F will measure \w) with probability 



Proof. The evolution of the Hamiltonian H is e - ltH / h = e *Et/h e -itH/h^ n 



4 Digitizing The GFG Algorithm 

In this section, we will digitize the GFG algorithm. The philosophy which 
we adapt here is that we are aware of the GFG algorithm and have the Lie- 
Trotter product formula at our disposal. We will see that under this philosophy, 
digitizing the GFG algorithm naturally leads to the Grover algorithm. For those 
who are familiar with Feynman path integral methods, it is also possible to view 
the Grover algorithm as a by product of giving the GFG algorithm a Feynman 
path integral interpretation. Of course this point of view is hindsight since the 
development of the two algorithms came in the reverse order and we are aware 
of Grover's algorithm. It is not surprising that this is hindsight since Feynman 
path integral methods are in general more complex than operator methods 4,1 . 
We take on the approach in this section because historically, quantum mechanics 
is a time dependent analog theory and because it is pedagogical in answering 
the question: what is a quantum algorithm? 

We start with the Lie- Trotter produce formula. For our purposes, we only 
need the finite dimensional version, which is the Lie product formula, see [20]. 

Theorem 4.1. Let A and B be self-adjoint, finite- dimensional matrices, then 



e -iv(A+B) = Um 

n^oo 

Furthermore, for any k, 



e -iT)A/n e -irjB/n 



(4.1) 



^ e - iv (A+B) _ [ e -i V A/k e -t n B/k^k^ < V 2 ° (II [A -B] II) ^ ^ 2 j 

where [A, B] = AB - BA. 

41 For example, it was many years after Fcynman's formulation of the Feynman path integral 
before the hydrogen atom path integral was solved (with physicist's rigor). Feynman himself 
was not able to solve the hydrogen atom path integral. See [13] and references within. 
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Proof. Let S k = e ~iv(A+B)/k^ and Tfe = e -v n A/k e -i n B/k _ Then 



1^-^11 = 



fc-1 



(4.3) 



^ST^-TfcJT*- 1 - 

m=0 

< fcll^rilll^- 1 - 1 "!!!! (5 fc - T k ) || = fc|| (S fc - T k )\\. 
Further, using the Campbell-Baker-Hausdorff approximation yields 

IK^-^II.S^M. ,4.4, 

Expressions 4.3 and 4.4 imply 

|| e - ir7 (A+S) _ [ e -i V A/k e -t n B/k^k^ = ||£fc _ TgW < — ^il^lllj _> o. □ 

(4.5) 

Theorem 4.1 is the first step towards a Feynman path integral interpretation 
for the GFG algorithm. We will derive a sum-over-paths interpretation in section 
5. We continue with our attempt to digitize the GHG algorithm. 

Proposition 4.1. Let P be a projection operator, then 

e- mP = I-2P (4.6) 
Proof. Power expanding and using the fact that P is a projection, we get 

(— i-wP) 71 (— i7r) n 



= 1 , ' = K —^P = I+{e-™-l)P = I-2P. □ (4.7) 



n— n— 1 



Proposition 4.2. Let = |tu)(io|, then Uf — e mPw . 

Proof. This is an application of proposition 4.1 by writing {// = / — 2P W . □ 
Proposition 4.3. Let P s = \s)(s\, then U s = e -W«--0. 

Proof. This is an application of proposition 4.1 by writing U s = —(I — 2P S ) = 

We are now ready to digitize the GHG algorithm presented in theorem 3.1. 
The Hamiltonian which we use for the GFG algorithm is H = E(P S — I + P w ), 
which produces the evolution e -«*- E (- p s-- f + p ™)/ s Using the Lie- Trotter product 
formula on this evolution yields 

e -itE(P e -I+P w )/h _ j im ( e -itE{P s -I)/nh e -itEP w /nh\ n _ ( 4 _g) 
n^oo \ / 

Further, for large k, we can approximate the evolution via 

e -itE(p s -i+p m )/h _ ^ e -itE(p s -i)/kh e -itEP a /kh^ k ^ 4g ^ 



8 



j.2 IT'S j.2 Tp 2 

where the error in the approximation is of order fe2 , . For k >> -p— , this 
approximation works well. Suppose we let ^ = fc7r, then according to propo- 
sitions 4.2 and 4.3, he right hand side of approximation 4.9 becomes (U s llf) k , 
which is a fc-product of the "Grover search engine". With t§- = kir, the error 
bound in approximation 4.9 is no longer in inverse powers of k, it is (1). This 
can be seen as follows. The Lie- Trotter product formula tells us that the error 
is bounded by 

t^O(\\[P s -I,P w }\\) 



k 

Further, 



(4.10) 



[P S -I,P W }\\ = 01 -L], (4.11) 



and for the GFG and the Grover algorithms, 

* = o (Vjv) , k = o (Vjv) . (4.12) 

Hence, 4.10 becomes (1). This result, which is not our goal, was obtained 
by [21]. Our goal here is to obtain error bounds in inverse powers of k. Thus, 
in the asymptotic limit on the size of the database to be searched, the success 
probability of the Grover algorithm approaches the success probability of the 
GFG algorithm. 

As far as running the GFG algorithm is concerned, we are interested in the 
probability 

T(t) = |H e -^| S >| 2 -sin 2 ^ +2/ 2 cos 2 ^, y = yi. (4.13) 

Thus, we settle for values of i's and fc's such that the following three properties 
are satisfied. First, t = (VN^J. Second, the right hand side of 4.9 becomes 

a fc-product of the "Grover search engine". Third, with those values of t's and 
fc's, the probabilities given by equations 4.13 and 2.8 are close to each other 
with error bounded by inverse powers of k. In other words, rather than globally 
digitize the GFG evolution operator, we settle for digitizing the operator in such 
a way so that the success probability of the GFG algorithm is well preserved. 
More simply, we wish to find correct combinations of t's and fc's so that 



1 w 3>(i) = \{w\e~ UH ' h \s)\ 2 = 



(w\ lim ( e -itE(P 3 -i)/nh e -itEP w /kn\ n \ s) 



(w\ (^ e - itE ( p s- r )/kh e -itEP m /kh^ k 



\{w\{U s U f f\s)\\ 

(4.14) 

where the error of the second approximation is bounded by inverse powers of k. 
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Theorem 4.2. Let i,meB! be such that l,m « VN, m is odd, and 

4tt/ 



7T - 2 



m + e, 



where \e\ ~ 0. Let k 
Finally, let 



2ttIVW 

tt-2 



be the nearest integer to 2rr J_^ , and let t 



Then, 



5 = max{e A l—y 



^ e -itE(p e -i)/kh e -itEP w /khy = (u s U f ) k , 



(4.15) 

~ E ■ 

(4.16) 
(4.17) 



p(t) = \(w\e- ttH/h \s}\ 2 



(4.18) 



and 



y k = 



,w e 



-itE(P,-I)/kh -itEP w /kh 



) l»> 



= \ (w\(U s U f ) k \s)\ 2 =cos 2 (kd + a) =y(t) + Q 



N 



(4.19) 



Proof. Equation 4.17 follows from propositions 4.2 and 4.3. Recall that (see [4]) 
'« I -= I = J + I VT? I . ( 4 - 20 ) 



= sm 



'N 

V N 

With the above values of fc, we have 



2 V^ + (^, ! 



Further, 
Et 



kn 2-rrlVN tt 



v^ft v/iV tt-2 V^V 
Then, 

M 4ttZ 2tt 2 Z 



+ Oh - 



ViVft 7T - 2 TT-2 



I TT 



-27T/ + OI \/ — I . 



(4.21) 



(4.22) 



(4.23) 



(4.24) 
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Hence, 



(to | (e 



itE(P s -I)/kh -itEP w /kh 



) l«> 



\(w\(U s Uf) \s)\ = cos 2 (k9 + a) 



cos 



= sin 



= 3>(t) + 




As for P(t), we have 




= sin 2 (| m + (5))+ o(l 



since to is odd, 4.26 becomes 



(4.25) 



sin2 f *M 



A' 



(4.26) 



P(t) = 1 + 0(6 2 ). □ 



(4.27) 



In light of theorem 4.2, we can interpret Grover's algorithm as a digitization 
of the GFG algorithm via the Lie product formula. While the wave function is 
not preserved under the digitization, the success probability is well preserved. 
We can also view the "Grover search engine" as e - ttE ( p s- -0 / ' kh e -*tEP w / 'kh w ;^j 1 
t and k as specified in theorem 4.2. Obviously, this view might only be useful 
for theoretical analysis since it would not be efficient to run Grover's algorithm 
by running t = (y^^j an d then k = (y^^j copies of the "Grover search 
engine" . 

We end this section with an numerical example of theorem 4.2. Suppose 
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iV = 2 30 , and/ = l. Then, 
4ttZ 



- 1 1 + . 00775358 = to + e, 

7T — 2 

S = max je, j = max |. 00775358, ^5 | = e, 
a = cos' 1 ( ^ ) w 1.57076581, 



/ 2\/2 30 — 1 \ 
6> = sin" 1 Q I w .00006104, 



2tt2 



15 



7T-2 

M fc?r 



180351, 
j 17.29093889. 



With the values given in 4.28, we have 



1 

.99985167 = 1 + (e 2 ) = 1 + (.00006012) . 



.99983048 = .99985167+0 (.00003052) = P(t) + 



N 



(4.28) 



P(t) w sin 2 (17.29093889) + ^ cos 2 (17.29093889) 



Lastly, 

P fc = cos 2 (fc0 + a) w cos 2 (180351 * .00006104 + 1.57076581) 

1 \ (4.30) 



Notice that in the example, we did not specify the value of E. Hence, we do 
not have a numerical value for t. 



5 Sum-Over-Paths Interpretation 

In this section, we will derive a sum-over-paths (Feynman path integral) in- 
terpretation for the GFG algorithm. Recently, sum-over-paths techniques were 
applied to quantum algorithms, see [6]. The sum-over-paths interpretation of 
the work in [6] is somewhat unconventional compared to traditional Feynman 
path integral methods (see [8, 9, 13, 22]) in the sense that the number of paths 
in [6] are finite. It is the goal of this section to derive a sum-over-paths interpre- 
tation for the GFG algorithm using more traditional methods. In traditional 
Feynman path integrals, breaking up the evolution operator e - rf ( H i+ H 2) v j a 
the Lie- Trotter product formula is the first step towards the derivation of the 
Feynman path integral. 

We are interested in deriving a sum-over-paths formula for the propagator 
of the GFG algorithm. For any quantum states \x) and \y), we can write 

{x\e- itH ' h \y) = Y,{Ame- im/h \k){k\y). (5.1) 
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The propagator K(j, k, t) is defined by 

K(j,k,t) = {j\e- itH/h \k). 



(5.2) 



The traditional non-relativistic quantum mechanics interpretation of the 
propagator K (x, y, t) is that it is the probability amplitude of a particle that 
starts at position x at time zero and ends up at position y at time t. The Feyn- 
man path integral interpretation of the propagator is that for any classical path 
from x to y at time zero to time t, an amplitude is assigned to that path and 
the total amplitude K (x, y, t) is the sum of all such path amplitudes. We are 
interested in obtaining a similar interpretation for the propagator in equation 
5.2. 



Proposition 5.1. Let \l) be a basis vector of the computational basis, then 

-UE\w)(w\ 



-itEf(l) 



Proof. We have 

-ME\w)(w\ 



io = io+E 



r=l 



itEY 1 
r 



n 



yP™|0 = \l)+(e=£r -l)5 Wll \l) 



+ {e=^ - l) f(l)] \l). 



(5.3) 



(5.4) 



Notice that 

{ —itE 

Hence the proposition follows. □ 
Proposition 5.2. Let \l), \m) be basis vectors of the computational basis, then 



- itEi ,}f,\ -I) -itE(-l) 

(l\e »« \m) = e ™ fi 



5 hm + ^(e-nn E -l) 



(5.6) 



Proof. Commuting the identity matrix and then Expanding the exponent yields 



,tEl <,{>\-I) -itS(-l) 

(l\e ^ |m) = e — — 



Sl,m + (e - l) (l\s)(s\m) 



its(-l) 



= e 



□ 



(5.7) 



Proposition 5.3. Let Iq = k and l n = j, then The propagator in equation 5.2is 
given by 



K(j,k,t)= lim £ [] 



-itE(f(l m )-l) 

e ™ ft 



(e^-l) 



il,; 2 |...,in-l ™=o 

-vn— 1 / ^/i \ -i\ n— 1 



lim 

n^oo 



E ex p 



Zi ,^2 ,...,i n -i 



-^E^(/('m)-l) 



n 

m=0 



TV 

(.*-.) 



(5.8) 
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Proof. Using the Lie- Trotter produce formula, equation 5.2 can be written as 



-itE(\3)(s\-I) -itE\»,)(w\ 

e " fi e 



K(j,k,t) = lim (j'l 
For each of the n products, inserting the identity 

^ = ^IU(U a = l,2,...n-l 



|fc). 



and denoting 
yields 

K(j,k,t)= lim 



n-1 



;i,i2,...,in-i m=0 



,+i|e 



-itE(| 3 )(a|-J) -itS| m )(™| 



(5.9) 



(5.10) 



(5.11) 



The result follows by applying propositions 5.1 and 5.2 to the last equation. □ 

Unfortunately, in the continuous limit, a sum-over-paths interpretation of 
proposition 5.3 is problematic. The expression <5; m ,; m+1 has the value 1 if at 
time mt/n, the state of the system is \l m ) and at time (m + l)t/n, the state 
of the system is |Z TO +i). In the limit to continuous time, 5i m j m+1 becomes 
meaningless. 

We proposed a nonstandard analysis formulation of proposition 5.3 which 
allows us to give proposition 5.3 a sum-over-paths interpretation. For details on 
nonstandard analysis, see [1, 5, 12, 23] and references within. 

Theorem 5.1. Let lo e *N be an infinite nonstandard natural number, then the 
propagator for the GFG algorithm is given by 



K(j,k,t) = st< 
J2 ex P 

Jl,Z2v)'w-l 



^Eri (/(M-i) 



n 

m=0 



(5.13) 



Proof. This is an application of nonstandard analysis's formulation of limits to 
proposition 5.3. □ 

The interpretation of theorem 5.1 is as follows. First, fix an infinite lo E *N 
and time slice [0, t] into lo + 1 number of steps T = {0, t/u, 2t/ui, . . . ,t}. Notice 
that e = t/u) is infinitesimal. For any path 



P:T-{1,2,...,JV} 



(5.14) 
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we associate the amplitude 



-^E^(/(fM)-i) 



W-l r 




exp 



m=0 L 



(5.15) 



to P. The propagator for GFG algorithm evolution is then infinitesimally close 
to the sum of the amplitudes of all possible paths P. Notice that for all fixed 
infinite uo e *N, the amplitude of any path has to be of the form given in 
equation 5.15. 

6 Semi-classical Limits 

In this last section, we will investigate the semi-classical limit behavior of the 
GFG and the Grover algorithms within the context of this paper and discuss 
some open problems. 

One way to obtain the semi-classical limit of a quantum system is by taking 
the limit ft — > 0. In the literature of quantum computing, ft is usually set to 
the value of 1. Throughout this paper, we have included ft in all the equations. 

According to section 4, the GFG algorithm is obtained by setting t = — — -. 

E 

Further, the two unitary operators Uf and U s were also obtained from the Lie- 
Trotter product approximation by setting t = ——. Hence, in the context of 

E 

this paper, in the semi-classical limit, both the GFG and the Grover algorithm 
must run instantaneously. This is because t — > as ft — > 0. It is not clear if this 
has any meaning in terms of the standard interpretation of Grover 's algorithm 
in the literature. In the literature, the unitary operators U s and Uf are time 
independent which could unrealistic since quantum mechanics is evolved by 
time dependent evolution operators. If we view decoherence due to interactions 
with the environment as the quantum system becoming classical, then we might 
consider modeling decoherence by ft — > 0. If we model decoherence this way, 
then the semi-classical limit behavior of the GFG and the Grover algorithm tells 
us that the algorithms must run instantaneously before the systems become 
decoherent. 

One could argue that semi-classical limit behaviors are not relevant to quan- 
tum computing, but recently there has been work done on the semi-classical 
limit behaviors of Shor's algorithm, see [10]. The authors in [10] proposed 
that it might be worthwhile to consider building quantum computers to run 
Shor's algorithm by using semi-classical devices. The authors showed that in 
the semi-classical limit, Shor's algorithm's success probability is not too severely 
modified. It would be interesting to see whether techniques similar to the ones 
used in this paper can be applied to Shor's algorithm and if in the semi-classical 
limit, Shor's algorithm must run instantaneously. 
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